Shear viscosity to entropy density ratio of a relativistic Hagedorn resonance gas 
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The new state of matter produced at Relativistic Heavy Ion Collider reveals a strongly coupled 
quark-gluon plasma with an extremely small shear viscosity to entropy density ratio rj/s. We 
calculate the rj/s of an equilibrated hadron matter characterized by a relativistic hadron resonance 
gas with a Hagedorn mass spectrum that grows exponentially with the hadron mass. We find with 
increase in temperature of the system the rj/s value decreases due to rapid increase in the multiplicity 
of massive resonances. In the vicinity of the critical temperature for deconfinement transition, the 
minimum value of rj/s in the Hagedorn resonance gas is found to be consistent with the current 
estimates for a strongly coupled quark-gluon plasma. 
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Heavy ion collisions at the BNL Relativistic Heavy Ion 
Collider (RHIC) have revealed a new state of matter [J- 
4] comprising of strongly interacting quarks and gluons 
(sQGP) This conclusion is based on viscous hydro- 
dynamic model analysis of elliptic flow that requires a 
very small shear viscosity to entropy density ratio rj/s 
of 0.08 — 0.24 [6-8]. The main uncertainty in these es- 
timates stems from the equation of state and the initial 
conditions employed. This value is remarkably similar to 
the lower bound rj/s > 1/4-7T obtained by Kovtun-Son- 
Starinets (KSS) .9] for infinitely coupled supersymmctric 
Yang-Mills gauge theory based on the AdS/CFT duality 
conjecture. A recent lattice calculation [lfj of gluonic 
plasma do support the current estimates. While leading 
order perturbative calculations result in a significantly 
large value of rj/s « 0.8 for a s = 0.3 (llj . 

It has been also argued 12, [l3| that with increasing 
temperature in the hadronic phase rj/s decreases and 
reaches a minimum at or near the critical temperature T c 
and increases thereafter in the deconfined phase. Indeed 
this behavior has been observed in several substances in 
nature all of which satisfy the KSS bound suggesting the 
bound could be universal. However, no consensus have 
yet been reached [l4[ on the physical mechanism that ac- 
counts for the thermodynamic and transport properties 
of the system that leads to the minimal viscosity. Under- 
standing the transport coefficients of the matter formed 
at RHIC is thus very challenging. Since in a heavy ion 
collision the system evolves from a QGP phase to the 
confined hadrons (low temperature phase of QCD), it 
is important to understand systematically the transport 
properties of a hadronic system in order to ascertain the 
sQGP properties with minimal uncertainty. In this letter, 
we investigate the transport properties of a infinite equi- 
librated hadronic matter comprising of massive hadronic 
resonances — the Hagedorn states 15-13] and all the low- 
lying observed hadrons within a Monte Carlo sequential 
binary emission model [3, [h| . 

Several calculations of rj/s in the hadron phase have 
been performed with various techniques such as the chi- 



ral per turbation theory, linearized Boltzmann equation 
[2CH221 ] . and microscopic transport theory 23]. However, 
most of these calculations employ at most two compo- 
nent system or mesonic gas with a wide variation in the 
estimate of rj/s ~ 0.08 — 1 at ~ T c . A recent microscopic 



transport calculation [24| within the UrQMD model ac- 
counted all the measured hadron species of mass to < 2 
GeV yields a minimum rj/s ~ 0.9 at zero baryon chemical 
potential. This is significantly higher than in the viscous 
hydrodynamic estimate. Clearly it indicates the impor- 
tance of massive resonances on the transport properties 
of hadron gas in the strong interaction domain of QCD 
matter especially near the critical temperature T c = 196 
MeV predicted in lattice simulations [251 ] . 

On the other hand, it was proposed [15j that the den- 
sity of hadronic states grows exponentially with reso- 
nance mass to, pw$(rn) ~ m~ a exp(m/Tjj), where Tg ~ 
150 — 200 MeV is the Hagedorn temperature. The mea- 
sured hadronic states up to m ~ 2 GeV are indeed qual- 
itatively consistent with the Hagedorn mass spectrum 
[261 ] . In absence of hadronic interaction, the energy den- 
sity for such a system will diverge leading to a maxi- 
mum/limiting temperature Tjj for the hadronic matter. 
In the infinite mass limit, the thermodynamical quan- 
tities for the Hagedorn resonance gas would thus show 
critical behavior as it crosses Th ~ T c which may be as- 
sociated with deconfining transition [27| • It has been also 
argued that the idea of Hagedorn states (HS) arise natu- 
rally in QCD at large N c [28[ and their indirect evidence 
can be found in lattice studies of thermodynamic proper- 
ties of hadron matter above the deconfinement temper- 
ature 29]. For a hadron resonance gas with Hagedorn 
states the reduction of rj/s in the vicinity of T c has been 
demonstrated [3fj| within the rate equation approach. 

In the present study of the properties of an equilibrated 
Hagedorn gas in the binary emission model [3, EH , we 
include the Hagedorn states (mass > 2 GeV) and all 
the low-lying measured hadrons available in the Particle 
Data Book. Detail description of the decay widths of HS 
and other resonances and their formation cross section in 
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FIG. 1: Time evolution of particle densities for an equilibrated 
system of hadron resonance gas with Hagedorn states at a zero 
net baryon density at temperatures T — 154 (top panel) and 
196 (bottom panel) MeV. 



binary collisions during_dynamic evolution of the system 



can be found in Ref. 



For brevity, we mention 
that depending on the mass, a HS may undergo binary 
decay into one of the three possible channels (i) two ob- 
served discrete hadrons (DH) of to < 2 GeV, (ii) a DH 
and a HS, (hi) two HSs. Based on Hagedorn hypothesis 
the density of massive states are assumed to be 

/ . , exp[{m - m g f(m - m g )} /T H ] 
p(m,q)=A- r J - ^ — , (1) 

{m — m g f(m — m g )} + m 2 , 

in the usual notation [l8],[__I_|. Here m r = 0.5 GeV, and a 
HS of mass to is characterized by its baryon, strangeness, 
spin and isospin quantum numbers q = (B, S, J, I, I z ) 
with a ground state mass m g = aQ(max|3i? + S\, 21) + 
as\S\; the parameters (a,Q,as) are determined empiri- 
cally from the measured smallest masses. This model was 
found to successfully describe the stable and resonance 
yield ratios and their spectra at RHIC. We consider here 
the Hagedorn temperature Tjj = 196 MeV consistent 
with the critical temperature T c in lattice QCD predic- 
tion of a crossover transition at vanishing baryon chemi- 
cal potential [___[. A resulting exponent a — 2.44 is then 
estimated by comparing the theoretical and experimental 
cumulants of the spectrum 18, 2__jj]. 

To generate the equilibrated matter we initiate a sys- 
tem of Hagedorn states in a box with periodic boundary 
conditions in configuration space. Subsequent binary de- 
cay and their regeneration in binary collision during time 
evolution enforce the system to thermodynamic equilib- 
rium. Chemical equilibrium is determined from the satu- 
ration of various particle densities when their production 
and annihilation rates become identical. While kinetic 



equilibration occurs when the system approaches momen- 
tum isotropization. By fitting the particle energy spectra 
with a Boltzmann distribution dNi/d 3 p ~ exp(— Ei/T*), 
thermal equilibration is verified when all the particle 
species at later times can be described by the same slope 
temperature T* . For point-particles employed in our sim- 
ulation, the pressure can be evaluated from the virial 
theorem as 



P P t(T*) 



1 



|P|: 



3V P 9 



(2) 



where and p\ are the momenta and energy of the ith 
particle. While the inclusion of Hagedorn states in the 
system provides an attractive interaction, use of point 
particles in the simulation ignores the repulsive interac- 
tions among the finite size hadrons especially the massive 
HS 17 1 that would drive the system to deconfinement 
at T c . The effects of repulsive interaction can be included 
via excluded- volume approach 16[ where the volume oc- 
cupied by ith hadron in the relativistic approximation 
is Vi = p®/4B, where B is the effective MIT bag con- 
stant. The excluded- volume temperature and pressure 
are related to their point-particle analogs as [l_i| 
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l-P P t(T*)/4B 



, P(T) 



P P t(T* 



l-P P t(T*)/4B 



(3) 



The other thermodynamic quantities can be obtained us- 
ing the thermodynamic identities. We consider B 1 / 4 — 

0. 40 GeV that results in P(T C ) < AB with a correspond- 
ing limiting temperature T lim > T c [l6[ as evident from 
Eq. ([3]). We restrict our calculations to T < T c in the 
confined hadronic system. 

Figure [T] shows the time evolution of various particle 
abundances at zero net baryon p_g and strangeness ps 
densities. At a temperature T w 154 MeV, the initial 
distribution of massive Hagedorn states (m ~ 20 GeV) 
decay and regenerate within few fm/c. Subsequently the 
HS multiplicity decreases rapidly to a small value with si- 
multaneous production of lighter (to < 2 GeV) hadrons. 
In contrast, at T ps 196 MeV each massive HS decays 
dominantly to a lighter HS and a discrete hadron result- 
ing in a continuous increase in the density of Hagedorn 
states and other hadron resonances. At about t = 200 
fm/c the particle densities saturate suggesting chemical 
equilibration has been achieved. 

In Fig. [2] (left panel) we present the equation of state, 

1. e. pressure versus energy density, of an equilibrated 
hadron matter at zero net baryon density. With in- 
creasing energy density when the temperature T ~ T c , 
enhanced particle production of especially heavy reso- 
nances and HS causes softening of the equation of state. 
Consequently the speed of sound c 2 s = (dP/de)v in the 
medium drops gradually from c 2 s pa 0.18 at T ps 100 MeV 
to about 0.09 at T c . The estimated minimum speed of 
sound is consistent with the lattice data 25|. In contrast, 
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FIG. 2: Pressure versus energy density (left panel), entropy 
density versus temperature (right panel), and the interaction 
measure Q = (e — 3P)/T 4 (inset). The results are in sequen- 
tial binary decay (SBD) model for hadron gas with Hagedorn 
states (HS) (circles) and in the statistical model (SM) for an 
pion gas (stars), a hadron resonance gas without HS (trian- 
gles), and with inclusion of HS (squares). 



the UrQMD model hadron resonan ce g as (without HS) 



24] 



yields a constant value of c? s = 0.18 

To gauge the dynamics of hadron-hadron interaction, 
we also investigate the thermodynamic properties in an 
independent statistical model for an hadron resonance 
gas in the grand canonical ensemble. The particle densi- 



(4) 



ties in the statistical model is given by [31 



47rp 2 dp 



2.3./ PM)dm^ 



exp(p? - m)/T* ± 1 



where the ± sign refers to fermions/bosons. For the 
ith species, gi is the spin-isospin degeneracy and fa — 
fj-sBi + fasSi is the chemical potential with ps and ps 
the baryon and strangeness chemical potentials, respec- 
tively. For the observed resonances, pi(m) is taken as 
Breit-Wigner mass distribution while for HS it is replaced 
by Phs °f Eq. (HJ). The results are found to be rather 
insensitive to the choice of upper mass limit of integra- 
tion for HS m max > 40 GeV. In the limit m max — > oo, as 
the estimated exponent a — 2.44 in Eq. ((T|), the point- 
particle energy density, pressure and entropy density of 
an ideal classical (Boltzmann) Hagedorn gas remain fi- 
nite at T < Th and diverge at T > Tjj exhibiting critical 
behavior [27J; albeit in the excluded- volume approach, 
the corresponding quantities evaluated at T < Th = T c 
would remain finite. 

In Fig. [2] (left panel) we also show the equation of state 
obtained in the statistical model. Compared to a pion gas 
(stars), inclusion of measured discrete resonances up to 
mass ~ 2 GeV (triangles) results in a considerable de- 
crease in the system pressure. The interaction measure, 



6 = (e - 3P)/T 4 , then gradually increases as T — > T c . 
Further inclusion of Hagedorn states in the statistical 
model opens up additional degrees of freedom and pro- 
vides an attractive interaction that slows the pressure 
increase at e > 1.5 GeV/fm 3 . The trace anomaly 8 then 
increases dramatically near T = Th ■ The speed of sound 
in the statistical model with HS turns out to be rather 



rail 



0.07 at T c . Comparison of this curve with se- 



quential binary emission model result provides a measure 
of the dynamics of repulsive interaction in the medium. 
The interaction effect starts at temperature above the 
pion mass and becomes significant at high energies where 
massive resonances are produced. 

The point-particle entropy density of the system can 
be computed using the Gibbs formula 

flpt (T*) = (e pt +Ppt-X>ft)/ T * 

i 

= (s P t + Ppt - pbpb)/T* , (5) 

where pg = BiPi is the net baryon density and the 
(initial) net strangeness density ps is set to zero in our 
calculation. The baryon chemical potential ps is evalu- 
ated from the particle yield ratios at equilibrium. Fig- 
ure [2] (right panel) shows the excluded- volume entropy 
density [16j as a function of temperature at zero baryon 
chemical potential. In the statistical model, compared to 
the pion gas, larger degrees of freedom in the hadron (and 
Hagedorn) resonance gas at the same energy density low- 
ers its temperature which results in higher entropy (32|. 
In the sequential emission model the increase (or slope) 
of entropy density with temperature near T c is somewhat 
less dramatic and captures the trend associated with that 
of a crossover transition obtained in the lattice data [2|| . 

We now extract the shear viscosity for the dynamically 
evolving system of an equilibrated hadron resonance gas 
with Hagedorn states using the Kubo relation [33[ ■ The 
transport coefficients determine the dynamics of fluctu- 
ations of dissipative fluxes about the equilibrated state. 
The Kubo formalism employs the linear response theory 
to relate the transport coefficients as correlations of dissi- 
pative fluxes by considering the dissipative fluxes as per- 
turbations to local thermodynamic equilibrium 33j-[35j. 
The Green-Kubo relation for shear viscosity is 



V = 



dt(ir x y(O,0)Tr x y(r,t)) cqlh . (6) 



Here t refers to time after the system equilibrates which 
is set at t = 0, and n xy is the shear component (traceless 
part) of the energy momentum tensor n^: 



d V —n-f(x,p) 



1 \ - PiPi 



V f-* 



=1 Pi 



(7) 



The equivalence accounts for the point particles used in 
our simulation. The averaging in Eq. ^ is over the 
ensemble of events generated in the simulation. 
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FIG. 3: Shear-stress correlation for the hadron resonance gas 
with Hagedorn states as a function of time at different system 
temperatures. The inset shows the temperature versus relax- 
ation time Ttt obtained from fits to the correlations using an 
exponential decay. 



Figure [3] displays the correlation function of the shear- 
stress tensor as a function of time. The correlation func- 
tions are found to damp exponentially with time, i.e. 
(71-^(0)71-^0)) - exp(-t/7v). Note at small tempera- 
tures (right panel) the correlations sustain for long times 
however show a power law tail. The relaxation times r w 
for the shear fluxes so obtained by fitting the correlation 
functions are found to decrease with increasing temper- 
ature. The shear viscosity is estimated from Eq. (JSJ by 
integrating the correlation function over time using the 
exponential decay found empirically. 
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FIG. 4: The shear viscosity to entropy density ratio rj/s as 
a function of temperature for a hadron resonance gas with 
Hagedorn states in thermodynamic equilibrium at various 
baryon chemical potential /j,b- The rj/s for a pure hadron res- 
onance gas in the UrQMD model [2J| is also shown (squares). 



The shear viscosity to entropy density ratio rj/s for the 
Hagedorn resonance gas in equilibrium at zero baryon 
chemical potential is presented in Fig. fU Also shown 
in the figure are the rj/s values for hadron resonance 
gas without Hagedorn states in the UrQMD model [24| . 
Our results for rj/s are comparable to the UrQMD esti- 
mates at temperatures T < 150 MeV. At higher temper- 
atures, massive resonances contribute dominantly result- 
ing in gradual reduction oirj/s with temperature [301 ] and 
reaches a minimum of (ry/s) m ; n w 0.34 at T c = 196 MeV. 
In contrast, the UrQMD model calculation shows almost 
a constant value for rj/s « 0.9 at T rj 125 - 166 MeV. 
The reduction of rj/s with Hagedorn mass spectrum can 
be realized as primarily due to enhancement of massive 
resonances that leads to a decrease in the mean free path 
of a particle and a corresponding increase in the average 
binary collision cross section a. In classical transport 
theory, as rj ~ p/a (where p is the average momentum of 
the particle), implies a reduction of shear viscosity. The 
(v/ s )mia ~ 0-34 at zero baryon chemical potential for 
the equilibrated Hagedorn resonance gas is significantly 
above the KSS bound of 1/4tt Q while somewhat close 
to the upper bound of 0.24 obtained from viscous hydro- 
dynamic analysis [7j, l8( of elliptic flow; albeit rj/s is as- 
sumed constant throughout the hydrodynamic evolution 
of the system. Note however, for Hagedorn resonance gas 
treated within the rate equation approach [30] , an upper 
limit of rj/s estimated in the kinetic theory is found to 
be as small as the KSS bound near T c . 

The measured antibaryon-to-baryon ratio B/B 0.8 
at midrapidity for central Au+Au collision at RHIC [3- 
4j suggests the hadronic matter formed possess, though 
small, but a finite net baryon number density ps- Dur- 
ing the late hadronic stage of evolution essentially from 
chemical equilibrium to thermal equilibrium the system 
could acquire non-equilibrium chemical composition [36j j . 
In fact, reproduction of the measured particle yield ratios 
at RHIC requires a non-zero baryon chemical potential 
in the sequential binary emission model [3 E3, while 



thermal model fits [37], |38( to the Au+Au collision data 
at the RHIC energy of y/s = 200 GeV yields a chemi- 
cal freeze-out temperature of T ~ 160 MeV and baryon 
chemical potential j±B ~ 20 MeV. In Fig. f3]we show rj/s 
as a function of temperature at a finite baryon chemi- 
cal potential [jis/T ss 0.13) created by inducing a finite 
Pb in the initial distribution of Hagedorn states. At a 
given T, the energy density of the system is larger for 
non-zero baryon chemical potential leading to further en- 
hancement of particularly massive particle abundances. 
In the T — jiB region explored, the shear viscosity rj in- 
creases compared to that at jis = 0. This increment 
can be understood classically — substantially large parti- 
cle densities though increase Oi and thereby reduces rji of 
each species i, however the total contribution, rj = rji, 
from these species is greatly enhanced [22j . Since the en- 
tropy density of the system also grows with T, the re- 



5 



suiting 7)/ s at ^b/T ~ 0.13 is found to be similar to 
/j,b = 0. In Fig. 2] we also illustrate the effect of high 
Hb/T « 3.0 associated with heavy ion collisions at the 
top AGS energy of ^/s = 4.85 GeV and at the lowest 



CERN/SPS energy of ^ = 6.27 GeV (371 |38|; albeit 
the temperatures reached in these reactions are much 
smaller at T ~ 125 - 145 MeV. We find that only for 
such large /is the extracted rj/s decreases further. The 
present study also suggests the need to incorporate, in 
the viscous hydrodynamics calculations, the Hagedorn 
states in the after-burning hadronic phase to quantify 
the dissipation and thus r\js ratio. 

While weakly coupled perturbative QCD estimate of 
i)/swlatr» T c , the QGP formed at RHIC in the re- 
gion T c < T < 2T C is thought to be strongly interacting. 

~ 3 



The current estimate of rj/s ~ 0.08 — 0.24 
at RHIC may be related to a minimum reached in the 
deconfined phase just above T c . Our estimate of rj/s for 
the Hagedorn resonance gas suggests that, as the sQGP 
cools, the rj/s in deconfined phase makes a relatively con- 
tinuous/smooth transition into the confined phase near 
T c without any discontinuity or a sharp increase proposed 
[241 ] in absence of Hagedorn states. 

In summary, we have studied the thermodynamic and 
transport properties of an infinite equilibrated matter 
composed of interacting hadron resonance gas plus an ex- 
ponentially increasing Hagedorn mass spectrum. At tem- 
peratures close to the Hagedorn temperature, Th =T C = 
196 MeV, the particle densities of the massive resonances 
are enhanced significantly. This leads to a sharp rise in 
the trace anomaly and the entropy density of the sys- 
tem which captures the trend seen in the lattice data at 
vanishing baryon chemical potential. The shear viscosity 
to entropy density ratio of the Hagedorn resonance gas, 
both in and out of chemical equilibrium, decreases with 
temperature and baryon chemical potential and reaches 
a minimum of (rj/s) m i n « 0.3 at about T c . The extracted 
minimum value is near the upper bound of current esti- 
mates for a strongly coupled quark-gluon plasma formed 
in ultra-relativistic heavy ion collisions at RHIC. 
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